NONTECHNICAL SUMMARY 


Resonant third-degree diurnal tides in the seas 
off Western Europe 


Does the ocean have normal modes? That is, does it have natural resonant 
oscillations like a vibrating bridge or a child on a swing? No ocean normal modes 
have ever been directly observed, but over the past two decades some indirect 
evidence for their existence has been slowly building up. The strongest evidence rests 
on observing a near-resonant response of the ocean to some large-scale forcing that is 
spatially well matched to a mode. 

This paper adds another piece of indirect evidence. In this case the large-scale 
forcing is a particular kind of tidal force called “third-degree diurnal.” Third-degree 
diurnal tides are extremely tiny, only a few millimeters or less, but in the seas off 
Western Europe they are seen to be significantly larger than expected (although still 
tiny). By studying this resonant amplification over a range of tidal frequencies, we are 
able to deduce the period and specific dissipation of the excited normal mode. 
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Abstract 


Third-degree diurnal tides are estimated from long time series of sea level 
measurements at three North Atlantic tide gauges. Although their amplitudes are 
only a few mm or less, their admittances are far larger than those of second-degree 
diurnal tides, just as Cartwright discovered for the Mi constituent. The tides 
are evidently resonantly enhanced owing to high spatial correlation between the 
third-degree spherical harmonic of the tidal potential and a near-diurnal oceanic 
normal mode that is most pronounced in the North Atlantic. By estimating 
the ocean tidal response across the diurnal band (5 tidal constituents plus nodal 
modulations), the period and Q of this mode and one nearby mode are estimated. 

1 Introduction 

This paper is a footnote to Cartwright’s studies of the curious Mj tide in the North At- 
lantic Ocean (Cartwright 1975, 1976; Cartwright et al. 1986). The Mi. constituent con- 
sists of a cluster of small spectral lines near 1 cycle per lunar day (period about 24.84 
hours). Table 1 lists the eight largest Mj lines in the astronomical tide-generating 
potential (Cartwright and Tayler 1971). The three largest lines are separated in 
frequency by only one cycle per 8.8 years. The middle of the three lines (155.555) 
is an exact subharmonic of the principal semidiurnal tide M 2 and arises from the 
third-degree term in the astronomical potential. Being dependent on the third, as 
opposed to the second, power of the moon’s parallax (the parallax is of order 1/60), 
third-degree tides are always very small, and in sea level data they are normally dom- 
inated by second-degree tides. Cartwright (1975) found the opposite to be true in 
the Northeast Atlantic where the ocean evidently magnifies the third-degree Mi line 
an order of magnitude relative to the second-degree lines (see Table 1). Cartwright’s 
explanation for this involves the ocean’s normal modes and their varying responses 
to different spatial forcing patterns. Since it forms the basis for the work presented 
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here, his explanation is worth recalling; it requires a brief digression into the theory 
of normal modes. 

A tide of frequency uj has a complex tidal admittance Z(uj, L) at location L which 
can be written as an expansion of oceanic normal modes (e.g., Platzman, 1991): 

Z(^L) = £(1 -u >/o k )~ l S k H k {L). (1) 

k 

According to this prescription each mode k contributes to the tidal admittance an 
amount depending on three factors: a “frequency factor” involving the mode’s com- 
plex frequency cq., a “shape factor” depending on the spatial coherence between the 
tidal potential and the mode elevation, and a location factor” which is simply the 
modal amplitude at L. . 

Normal modes for the world ocean have been computed by several groups. Ac- 
cording to calculations by Platzman et al. (1981) there are three modes within or very 
near the diurnal band, having periods 23.7, 25.7, and 28.7 hours. The latter mode 
is primarily an Antarctic Kelvin wave and Pacific half-wave, with hardly any influ- 
ence in the Atlantic Ocean In contrast, the other two modes are primarily Atlantic 
modes. The relevant portion of the 25.7-h mode in shown in Figure 1. It has two 
amphidromes rotating cyclonically in each basin and largest amplitudes in the North 
Atlantic, especially the northeast sector. The 23.7-h mode is not unlike the 25.7-h 
mode in being largest in the North Atlantic and relatively small elsewhere; it does, 
however, have an additional amphidrome in the equatorial Atlantic. See Platzman et 
al. (1981) and Platzman (1985) for a comprehensive set of mode diagrams. 

As Cartwright (1975) pointed out, the 25.7-h mode is spatially well correlated with 
the third-degree component of the diurnal tidal potential, but not with the second- 
degree component. The relevant spherical harmonics of the potential are (Cartwright 
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and Tayler, 1971): 


Y 2 (0, <f>) = y/(5/2‘in)3 sinO cos6 e 1 ^ 

= ^(7/192^)3 sin <9(5 cos 2 9- l)e i0 

for spherical coordinates (0,0). Both functions are of order 1, as is appropriate for 
diurnal tides; they are shown in Figure 2. The Y3 harmonic, like the 25.7-h mode, 
is symmetric about the equator, and both the harmonic and the mode have local 
maxima near latitudes 0° and ±60°. In contrast, the harmonic is anti-symmetric 
about the equator, and it cannot be expected to easily excite a symmetric mode. 
Platzman (1984) confirmed this argument by calculating the correlations: he found 
the shape factor S* is 0.30 for Y% but only 0.05 for Y^- (Note that the faster 23.7-h 
mode is also symmetric, but with an additional amphidrome near the equator it is 
less well matched to Y^ 1 than is the 25.7-h mode.) 

If this story of Mi is correct, then one naturally suspects that other third-degree 
diurnal tides are similarly enhanced in the Northeast Atlantic. It is the purpose of 
this note to show that they indeed are. In addition, by studying the third-degree 
tidal admittances across the diurnal band, one may deduce further characteristics of 
the normal modes and their responses, namely their periods — to be compared with 
Platzman’s — and their specific dissipation Q ~ l . 

Like Mi the other third-degree diurnal tides are separated from second-degree 
tides by only 1 cycle in 8.8 years. The other third-degree tides are also significantly 
smaller than Mi, which is itself only a few millimeters. Both facts suggest that very 
long time series of hourly data are required before reliable estimates can be obtained. 
The data employed below, from three tide gauges, span at least two complete cycles 
of the lunar node, i.e. at least 37 years. 
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2 Estimates of degree-3 tides 


Hourly sea level data were analyzed from three stations: Newlyn, England (50° 6' N, 
5°33 / W), Lerwick, Scotland (60°6 ; N, 1°8'W), and Vigo, Spain (42°14'N, 8°44'W). The 
data timespans are: 1960-1997 (Newlyn), 1959-1997 (Lerwick), 1943-1990 (Vigo). 
The data are of high quality and remarkably complete; the longest gaps are a three- 
week period during 1977 at Lerwick and two two-week periods during 1953 at Vigo. 

The tidal analyses were performed by least-squares fits to sinusoids of known fre- 
quencies. All major constituents, including long-period tides and nonlinear compound 
tides, were accounted for, although in the tables given here (Tables 2-4) only the rel- 
evant third-degree diurnals are listed. For these third-degree diurnals, solutions were 
attempted for all lines whosp potential amplitudes exceed 0.3 mm in the Cartwright- 
Tayler expansion. There are five such constituents spanning the diurnal band between 
2Qi and Ji. 2Qi is more than 10 times smaller than Mi, but it is included because 
of its importance at the edge of the band. The five constituents are separated in 
frequency by one cycle/month, with a single gap where one normally expects the Ki 
constituent (its largest third-degree line is only 0.22 mm, even smaller than 2Qi). 

Because of the extraordinarily long lengths of the time series, relevant nodal side 
lines in constituents can be included as separate estimates, independent of the main 
lines. For the third-degree diurnals, nodal lines above the 0.3 mm cutoff exist for 
constituents Qi, Mi, and J|. The nodal lines of Qi and Ji appear to yield genuine 
third-degree estimates, but the nodal lines of Mi are problematic: one of them clearly 
cannot be separated from second-degree lines, since (see Table 1) lines 155.565 and 
155.645 are of almost identical frequencies, and the other line may be affected by 
nodal modulations in a weak nonlinear interaction between L 2 and Ki (equivalent to 
a line at 155.465). Only nodal lines from Qi and Ji are included in Tables 2-4. 

The tables give amplitude and Greenwich phase lag and the magnitude of the 
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admittance Z, defined in this context as the ratio of the (complex) elevation to the 
normalized potential. All three quantities include standard error estimates, based on 
the diagonal elements of the least-squares covariance matrix. The admittance errors 
for 2Q! are understandably largest owing to its small amplitude in the astronomical 
potential. The given error estimates appear realistic according to two tests: (1) 
The original time series were partitioned into year-long segments and the statistical 
variability of the yearly estimates was compared with the tabulated standard errors 
(allowing for a >/n factor for n yearly estimates); this test, of course, can only be 
performed for major second-degree constituents, but the agreement in these terms 
appeared satisfactory. (2) The error estimates appear to account adequately for the 
small differences in \Z\ and Greenwich phase lag between the Qi and Ji nodal lines 
and their main lines, in keeping with a smooth oceanic admittance. 

The tidal estimates for Mi at Newlyn and Lerwick agree well with those obtained 
by Cartwright (1975) for these same two stations, and station Vigo similarly shows 
an enhanced Mi admittance. In fact, all three stations display a strong enhancement 
across the entire diurnal band. This lends considerable weight to Cartwright’s original 
explanation of a resonance between the tidal potential and a single operating normal 
mode (of Q not too high). Moreover, all five constituents show phase lags increasing 
towards the north, consistent with the normal mode phases seen in Figure 1. Both 
the enhanced amplitudes and the consistent phases must be counted as yet further 
evidence for the existence of oceanic normal modes, supplementing various other 
evidence discussed by Platzman (1991). 

At all three stations, the admittance at Ji is comparable to or slightly larger 
than that at M\. More interestingly, all three stations show admittance amplitudes 
increasing as the frequency decreases away from Mi, with 2Qx having the largest 
admittance in all cases. If Platzman’s single 25.7-h mode is the primary contributor 
to these tides, which is the case according to his tidal synthesis (Platzman 1984, Figure 
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10), then the admittances strongly suggest that his period for this mode is too low. 
A period below that of Qi, and possibly below that of 2Qi, appears more consistent 
with the observed tidal admittances. An error of 10-20% in Platzman’s period is not 
unreasonable given various approximations made in the original calculations (e.g., 
neglect of self-attraction and crustal loading). 

3 Modal fits 

With the tidal admittance defined across the diurnal band, it is straightforward to 
use the formalism of Equation (1) to fit hypothetical normal modes (e.g., Platzman 
1991). In particular, if one dominant normal mode is presumed to account for the 
behavior across the band, then we may fit the admittance Z by a simple equation of 
form 

Z(uj) = (1 — w/cr) -1 (2) 

where the modal strength b and frequency a are treated as free parameters. Here 
Z, b, and a are complex, while oj is real. Of primary interest is the frequency a 
which determines the modal period T = 27r/Re<7 and specific dissipation Q -1 = 

2 Imtr/Re<7. 

Figure 3 shows the result of fitting the tidal admittances at Newlyn, using a 
weighted, nonlinear, least-squares algorithm. In this case, there are 2 complex pa- 
rameters to be determined from 5 complex admittances. The resulting frequency 
corresponds to a modal period of 28.1 hours and Q of 7. This estimate of Q is about 
twice the value determined for the diurnal-band Antarctic Kelvin wave but less than 
the Q determined for the semidiurnal tides of the, North Atlantic (e.g., Platzman 
1991). It is clear that the fitted period depends sensitively on the admittance at 2Qi; 
since its amplitude falls well short of the linear trend between Oi and Qi, the fitted 
curve is forced downward, causing a peak (and hence estimated period) in the general 
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vicinity of 2Q[. In this region of the diagram, the fit is satisfactory given the sizes of 
the error bars. On the other hand, the amplitude fit at Ji is clearly poor and suggests 
the presence of one or more additional modes. 

Equation (2) has therefore been augmented with a second oscillator. The resulting 
fit to the diurnal admittances is shown in Figure 4. The estimated periods in the case 
are 27.5 and 23.0 hours, somewhat closer to Platzman’s 25.7 and 23.7 hours. The 
corresponding Q, around 10, is higher than the single-mode fit. One should not, 
however, put too much credence to this augmented fit; we are fitting 4 complex 
parameters to only 5 complex admittances, and the very good agreement shown in 
Figure 4 in fact suggests overfitting. Unfortunately, there are no additional tidal 
frequencies that can be used to constrain a two-mode fit. 

t 

Fits to the data at Lerwick and Vigo are quite similar to those at Newlyn. The 
estimated periods and Q values for all stations are summarized in Table 5. I am 
disinclined to make a choice as to the “best” estimated period for the fundamental 
mode, except to say that the available evidence clearly suggests a period longer than 
Platzman’s original 25.7 hours, a period most likely longer than the period of Qi, and 
perhaps in the vicinity of 27 to 29 hours. 

Acknowledgments. The tide gauge data employed here were obtained courtesy of 
the British Oceanographic Data Centre, Bidston Observatory, and the University of 
Hawaii Sea Level Center. 
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TABLES 


Table 1: Mj tidal amplitudes at Newlyn, England 


Doodson 

number 

Degree 

Frequency 

(°/h) 

Potential* 

(mm) 

Equilibrium** 

(mm) 

Observed*** 

(mm) 

155.445 

2 

14.48520 

1.37 

0.359 


155.455 

2 

14.48741 

7.41 

1.943 


155.545 

3 

14.48985 

0.59 

0.188 


155.555 

3 

14.49205 

3.99 

1.268 

4.8 

155.565 

3 

14.49426 

0.52 

0.165 


155.645 

2 

14.49449 

0.59 

0.155 


155.655 

2 

14.49669 

20.62 

5.407 

1.2 

155.665 

2 

14.49890 

4.14 

1.086 



Potential normalized by acceleration of gravity, extracted from tables of Cartwright 
and Tayler (1971). 

** Equilibrium amplitude equals potential x y n where y n accounts for elastic body 
tide. 

*** Sea level amplitudes estimated by Cartwright (1975). 
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Table 2: Third-degree diurnal tides at Newlyn. England 


Tide 

Doodson 

number 

Frequency 

(°/h) 

Potential* 

(mm) 

Amplitude 

(mm) 

Phase lag 
(deg) 

\z\ 

2Qi 

125.655 

12.84964 

0.35 

0.81 ±0.22 

205 ± 16 

2.31 ±0.62 


135.545 

13.39181 

0.50 

1.00 ±0.22 

247 ±12 

1.99 ±0.43 

Qi 

135.555 

13.39402 

1.28 

2.82 ±0.22 

249 ±5 

2.20 ±0.18 

o x 

145.655 

13.94768 

0.65 

0.93 ±0.22 

276 ± 14 

1.42 ±0.33 

Mi 

155.555 

14.49205 

3.99 

4.67 ±0.22 

275 ±3 

1.17 ±0.06 

Ji 

175.555 

15.59008 

1.46 

1.84 ±0.22 

294 ±7 

1.26 ±0.15 


175.565 

15.59229 

0.59 

0.77 ±0.22 

325 ± 16 

1.30 ±0.37 


* Potential normalized by acceleration of gravity, extracted from tables of Cartwright 
and Tayler (1971). 


Table 3: Thirty-degree diurnal tides at Lerwick, Scotland 


Tide 

Doodson 

number 

Frequency 

(°/h) 

Potential 

(mm) 

Amplitude 

(mm) 

Phase lag 
(deg) 

\z\ 

2Qi 

125.655 

12.84964 

0.35 

1.35 ±0.18 

266 ±8 

3.84 ±0.51 


135.545 

13.39181 

0.50 

1.66 ±0.18 

307 ±6 

3.31 ±0.35 

Qi 

135.555 

13.39402 

1.28 

4.08 ±0.18 

307 ±3 

3.18 ±0.14 

Oi 

145.655 

13.94768 

0.65 

1.66 ±0.18 

337 ±6 

2.54 ±0.27 

Mt 

155.555 

14.49205 

3.99 

7.25 ±0.18 

332 ±2 

1.82 ±0.05 

Ji 

175.555 

15.59008 

1.46 

2.54 ±0.18 

354 ±4 

1.74 ±0.12 


175.565 

15.59229 

0.59 

1.22 ±0.18 

15 ±9 

2.06 ±0.30 


Table 4: Third-degree diurnal tides at Vigo, Spain 


Tide 

Doodson 

number 

Frequency 

(°/h) 

Potential 

(mm) 

Amplitude 

(mm) 

Phase lag 
(deg) 

\z\ 

2Qi 

125.655 

12.84964 

0.35 

1.08 ±0.15 

179 ±8 

3.09 ±0.42 


135.545 

13.39181 

0.50 

1.02 ±0.15 

229 ±9 

2.04 ±0.30 

Qi 

135.555 

13.39402 

1.28 

3.11 ± 0.15 

213 ±3 

2.43 ±0.12 

O! 

145.655 

13.94768 

0.65 

1.29 ±0.15 

236 ±7 

1.98 ±0.23 

M! 

155.555 

14.49205 

3.99 

5.48 ± 0.15 

232 ±2 

1.37 ±0.04 

Ji 

175.555 

15.59008 

1.46 

1.83 ±0.15 

250 ±5 

1.25 ±0.11 


175.565 

15.59229 

0.59 

0.68 ±0.15 

264 ± 13 

1.15 ±0.26 
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Table 5: Period and Q from modal fits 


Station 

One-mode fits 


Two-mode fits 


Period 

Q 

Period 

Q 

Period 

Q 

Newlyn 

28.1 

6.9 

27.5 

11.4 

23.0' 

9.2 

Lerwick 

28.2 

6.3 

27.6 

9.2 

23.4 

9.9 

Vigo 

28.7 

6.7 

27.8 

6.9 

24.1 

11.2 


Periods are in hours. 
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FIGURE CAPTIONS 


Figure 1: The 25. 7- hour mode in the Atlantic Ocean as calculated by Platzman et 
al. (1981). Solid lines are 60° phase contours with arrowhead indicating direction 
of rotation. Dashed lines are amplitude contours, with shading denoting areas with 
amplitudes exceeding the global rms amplitude. Reproduced from Platzman (1991) 
with permission. 


Figure 2. Contours of the second-degree (top) and third-degree (bottom) spherical 
harmonic components of the diurnal tidal potential. Solid and dotted lines denote 
regions of opposite sign. 


Figure 3: Amplitude (top) and phase (bottom) of the third-degree diurnal tidal ad- 
mittance at Newlyn, England, corresponding to constituent estimates tabulated in 
Table 1. The smooth curve represents a one-oscillator fit to the data. The amplitude 
peaks at the modal period 28.1 hours. The poor fit to Ji presumably indicates the 
influence of one or more additional modes. 


Figure 4: As in Figure 3 but for a two-oscillator fit to the data. 




Figure 1: The 25.7-hour mode in the Atlantic Ocean as calculated by Platzman et 
al. (1981). Solid lines are 60° phase contours with arrowhead indicating direction 
of rotation. Dashed lines are amplitude contours, with shading denoting areas with 
amplitudes exceeding the global nna amplitude. Reproduced from Platzman (1991) 
with permission. 
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